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ABSTRACT: We present the first complete calculation performed within the Four Dimen- 
sional Regularization scheme (FDR), namely the loop-induced on-shell amplitude for the 
Higgs boson decay into two photons in an arbitrary gauge. FDR is a new technique 
-free of infinities- for addressing multi-loop calculus, which automatically preserves gauge 
invariance, allowing for a 4-dimensional computation at the same time. We obtained the 
same result as that assessed in dimensional regularization, thereby explicitly verifying, in a 
realistic case, that FDR respects gauge invariance. 
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1 Introduction 

Calculating higher order corrections in Quantum Field Theories (QFTs) is becoming more 
and more influential both from a phenomenological point of view and from a more theoretical 
perspective. The success of the Standard Model (SM) in describing the observed high energy 
data, especially in the light of the discovery of a Higgs boson candidate [1, 2], leaves little 
room for large and unexpected signals at colliders, at least at the present energy regime. 
One is therefore forced to rely on the full predictive power of QFTs and compute the tiny 
effects induced by the Radiative Corrections (RCs), in the hope of finding small deviations 
from the observed data. 

The computation of RCs is very demanding from a technical point of view. Methods 
used in the past decades for processes with moderate multiplicities are becoming less and 
less adequate in view of the growing complexity of the interesting final states. While a 
lot of work has been recently devoted to deal, in an efficient way, with 1-loop processes at 
large multiplicities [3-6] , very little simplification has been achieved so far in the field of the 
multi-loop calculations 1 . What seems to trigger most of the difficulties is the usual treat- 
ment, in the framework of Dimensional Regularization (DR) [11], of the infinities arising in 
the intermediate steps of the calculation: DR forces a huge and cumbersome analytic work 
to be done to ensure all needed cancellations, even before starting the actual computation of 
the relevant physical (finite) contribution. Therefore it looks obvious that major simplifica- 
tions are expected when performing the computations directly in the physical 4-dimensional 
Minkowsky space, also in the presence of Ultraviolet (UV) and Infrared/Collinear diver- 
gences. This observation led to several attempts to find four-dimensional alternatives to 

1 See, however, [7-10]. 
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the DR treatment of the UV infinities, such as differential renormalization [12], constrained 
differential renormalization [13, 14], which both work in the coordinate space, implicit renor- 
malization [15, 16] and LR [17], directly applicable in the momentum space. In a recent 
work, the FDR approach [18] has been proposed by one of us in which the UV problem is 
solved by simply re-interpreting the loop integrals appearing in the calculation. They are 
defined in such a way that infinities do not occur, at the price of introducing an arbitrary 
scale, called /i, which plays the role of the renormalization scale. The procedure works 
because the FDR re-interpretation respects gauge invariance by construction. 

The purpose of this work is to illustrate some of the assets of the FDR scheme, and to 
show how to put it into practice. To this goal, we present the first application of FDR to 
a complete calculation in the EW theory. In order to keep the most general approach, we 
have decided to work in an arbitrary i?£-gauge, thereby explicitly verifying that the method 
respects gauge invariance. We have chosen to compute the 1-loop on-shell amplitude for 
the Higgs boson decaying into two photons, which is known since a long time [19-22] and, 
given its relevance and simplicity, has been recently reconsidered by several authors as 
a case study [23-33]. Because there is no I/77 interaction at tree level in the SM, the 
process is finite; although no renormalization is needed, infinities arise at intermediate 
steps of the calculation, which demands to work within a divergence-safe framework, such 
as FDR. Moreover, being indirect, the process is mediated by either bosonic or fermionic 
loops. Calculating the bosonic contribution stands as a strong test of the gauge invariance 
property of FDR. On the other hand, the fermionic contribution gives us the opportunity 
of illustrating FDR in the presence of fermionic loops. 

The outline of the paper is as follows: in section 2 we review the FDR method, in 
section 3 we present tensor reduction in FDR and, in section 4, we report on the calculation 
details relevant for the process at hand. 

2 The FDR method in a nutshell 

FDR is a new method -gauge-invariant, free of infinities and 4-dimensional- to tackle loop 
calculations in QFTs. Here we review its main features and definitions, and refer to [18] for 
a more detailed discussion. 

Starting from the observation that UV divergences are unphysical and that they can 
be decoupled from the physically-relevant information of an amplitude, a new and coherent 
definition of loop integrals is possible, such that infinities are removed straight away and 
the integral is finite in 4-dimensions, with no need to absorb divergent quantities into 
renormalized parameters. Consider, for example, the dimensionally regulated integral 





where n = 4 + e is the number of space-time dimensions and 




with po = . 



(2.2) 



By introducing a new arbitrary scale /i, and defining 

Di = Di- fi 2 , q 2 = q 2 - fj 2 , 



(2.3) 
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the logarithmically divergent part can be separated from the the physical one with the help 
of the following partial fraction identity 



yielding 



lim 



d n q 1 



1 



lim 

At-5-0 




(2.4) 



d n q 



+ 



+ 



d 



lfD 1 q 2 Dq Di 



(2.5) 



where the divergent piece is written between squared brackets. The /u scale is essential 
in order to regularize the integral in the low energy regime, and it is understood as the 
scale with respect to which loop momenta are very large: this is why infinities are expected 
to decouple in the limit [i — > 0. Notice that the divergent contribution, being process- 
independent, effectively behaves like a vacuum bubble, and as such should not be taken 
into account when calculating a physical observable. The FDR integral corresponding to 
the integrand of eq. (2.1) is then defined to be 



? FDR 



= lim 



di 



dn 



d 4 g — ^+ ~"_ 



(2.6) 



so that it only includes the contribution conveying the physical information. The final iden- 
tification [i = fj,R effectively eliminates the dependence on the original cut-off, as explained 
in [18]. This reasoning can be generalized to I loops. Starting from a generic DR regularized 
integral 



^DR 



»R £ [fld n qi J({q 2 }), 

J i=l 



(2.7) 



where {q 2 } denotes the set of the squares of all possible combinations of loop-integration 
variables, the integrand J = Jy + Jp can be parametrized in terms of [i and split into 
vacuum configurations (Jy) an d finite part (Jf), so that one can write the FDR integral 
over J as 



a FDR 



i i 
[ Hld%] J({q 2 }) = lim f H d% J F ({q 2 }) 
J i=l M J i=l 



(2.8) 



where the replacement 

{q 2 }^{q 2 } (2.9) 

is performed in both numerator and denominators. The compact notation [ d 4 gj] in eq. (2.8) 
therefore implies: 

1. parametrizing in terms of 

2. separating and subsequently dropping the vacuum configurations; 

3. integrating in 4 dimensions; 



-3- 



4. performing the limit fM — y until a logarithmic divergence is met; 

5. evaluating the result in fi = jjlr. 

Note that one can also define the FDR integral as the difference between the DR integral 
and its vacuum configurations: 



from which it is manifest the invariance under the shift of any integration variable. 

Let us say a few more words on gauge-invariance, which is one of the key features 
of FDR, distinguishing it from other 4-dimensional methods. While it is obvious that 
terms conveying the kinematical dependence are equivalent in FDR and DR, a potential 
ambiguity remains in the constant term, because FDR subtracts vacuum configurations 
before integrating, while DR takes away poles in e after the integration. In general the 
constant term can be fixed by enforcing gauge invariance as an extra constraint of the 
amplitude. However, FDR and DR alike automatically respect gauge invariance, thereby 
leading straightforwardly to the same correct constant. In order for gauge-invariance to 
be preserved in FDR, it is crucial that the parametrization in terms of fj, is performed 
correctly, that is propagators Di and all q 2 appearing in the numerator should be promoted 
to their barred counterpart. This prescription is referred to as global treatment. A q and its 
associated /i should never be treated separately, which allows for the usual simplifications 
between numerator and denominator to take place. There are cases in which a fi 2 does 
appear alone (we will discuss an example in the next section). In this type of terms, fj? 
effectively plays the role of a squared loop momentum, and as such it should be treated. 
For example, the integral 



should be regarded as a logarithmically divergent one, that is to say it requires the same 
expansion needed to subtract the vacuum bubbles from 




(2.10) 
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namely 
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giving 



a FDI V) = iim f / A ( 2 - do - - + + r^=- 

M-xT J H \q 2 D D 1 D 2 q D\ D 2 q 6 D 2 



) 



(2.14) 
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The last integral behaves as l/// 2 , so that, in the limit \i — > it does not vanish, but 
generates a finite constant: 

3 FDI V) = ^. (2-15) 

This is equivalent to DR, when a finite term is obtained as the product of an 0(e)-term 
and a single pole. More explicitly one can prove that 

where q 2 = (q( n ') 2 — q 2 is the e-dimensional part of an n- vector, and the superscript (n) 
denotes an object living in n dimensions. 

As far as strings of Dirac matrices are concerned, the replacement rule q 2 — > q 2 is 
meant to be performed after calculating the trace (in 4 dimensions). The same result can 
be obtained by promoting ^ — > d = q ± // directly in the string 2 . This is achieved by 
defining ^ according to its position: 

(. . . j r 1 ■■ i an j ...) = (••• a ± ») t 1 • • • 7 Q -(^ t (-)» . . .) • (2.i7) 

The sign within the first ^ is chosen arbitrarily; in the following ^, the sign is opposite if an 
even number of 7-matrices occur between the two ^'s, and it is the same in the case of an 
odd number of 7-matrices. Likewise should be treated the subsequent pairs of ^'s occurring 
in the string. If chirality matrices are also involved, a gauge invariant treatment [34] requires 
their anticommutation at the beginning (or the end) of open strings before replacing ^ — > d. 
In the case of closed loops, 75 should be put next to the vertex corresponding to a potential 
non-conserved current. This reproduces the correct coefficient of the triangular anomaly, 
as observed in [18]. 



3 Passarino-Veltman reduction in FDR 

The main asset of the Passarino-Veltman [35] method (PV) is that it provides a gauge- 
invariant decomposition of the amplitude: after the reduction, the full amplitude is ex- 
pressed in terms of scalar integrals depending on physical thresholds only -multiplied by 
gauge independent coefficients- plus a rational part. 

In this section we review the PV reduction method within FDR. Our main aim is to 
illustrate the source of the constant terms contributing to the amplitude. For simplicity, 
we explicitly work out the decomposition of a rank-2 bubble. In the very same way higher 
rank tensors with more denominators can be decomposed, although, for brevity, we do not 
report here the result. 

The integral 

B^= [ [d 4 g]Xt, (3.1) 
2 Thanks also to the fact that FDR integrals involving odd powers of fi vanish [18]. 
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can be decomposed in terms of rank-2 symmetric tensors depending on p = p\\ 

B^ = B 00 g^ + B llP ^. (3.2) 

By contracting both sides with and p^Pv, we obtain a system of equations from which 
the scalar coefficients Bqq and B\\ can be extracted: 



g^ B^ = 4B 00 +p 2 B u 
P»Pv =p 2 (B 00 +p 2 B u ). 



(3.3) 



With the usual algebraic manipulations we can express the left hand sides of eq. (3.3) 
in terms of scalar integrals with 1 or 2 internal legs. Special care should be paid when 
contracting B^ u with the metric tensor; indeed, 



[d 4 d^W 



D Q D 1 



D D 1 

q 2 takes part in the usual simplifications between numerator and denominator, 
g, u = A (M 2 ) + M 2 B (p 2 ; M 2 , M 2 ) + / [d 4 q]^=- , 



(3.4) 



(3.5) 



while a fjr is left out. The last integral should be treated as that in eq. (2.11) by expanding 
— 1 — to extract all the terms behaving like i. Alternatively, one can use the definition of 
the FDR integral as a difference of integrals, as in eq. (2.10) 3 : 



lim / d 4 q < ^ 

J [ D D\ 



where 



D D 1 
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D D 1 
do + di 



+ 4 



{q-pf 



yielding 



[d 4 g]=L 
DoDi 



lim Id q n 



4„ ,.2 j dp + di | ^q-p) 2 



By solving the system in eq. (3.3), we finally obtain 
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+ 



-A (M 2 ) [p 2 - A_] + MM?) [2 P 2 - A. 



3p 4 



+ 



B (p 2 ;MlM 2 )[(p 2 -A_) 2 -p 2 M l 
3p 4 



(3.6) 



(3.7) 



.(3.8) 



(3.9) 



3 The integration can be performed in 4 dimensions because the divergent integrands cancel in the 



difference. 
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ki, n 



k 2 , v 



Figure 1. Sketch of the Feynman diagrams contributing to the process; the loop can be either 
fermionic or bosonic (see figures 2 and 3). The momenta are considered to be all incoming, and the 
virtual loop momentum q is outgoing from the Higgs vertex. 

where 

A ± = M 2 ±Af 2 . (3.10) 

This result is consistent with that obtained with the standard PV reduction in DR. 

Finally we point out that reduction methods other than PV can be used as well, e.g. 
the OPP method of [3]. The basic observation is that any algebraic manipulation of the 
integrand is legal in FDR, the only subtlety being the replacement of eq. (2.9) in the 
numerator, which may generate a rational part 4 . However, thanks to the correspondence 
in eq. (2.16), this contribution can be reinserted back -in an OPP like reduction of the FDR 
integrals- by using the same set of effective Feynman rules computed in DR [36, 37], or with 
the technique described in [38]. 

4 The Higgs boson decay into two photons 

Let us now concentrate on H — > 77. As mentioned above, because there is no #77 vertex 
in the SM Lagrangian, this process is finite and loop-induced. The generic amplitude is 
given by a bosonic and a fermionic contribution, which are independent of each other and 
separately gauge-invariant: 

M(M = (m w {P) + Y. n &) M f(v)) , (4.1) 

/ 

where /3 and i] are dimensionless kinematic parameters defined as 

4 M?,, 4 m 2 f 

By denoting with k\ and k 2 the momenta of the photons, as in figure 1, the amplitude 
reads 

M= M^E^eKh). (4.3) 

The tensorial structure of Ai^ u is dictated by on-shellness and gauge invariance, i.e. 
ki ■ e(ki) = and k*f M^ v = k^ M^ v = 0, so that 

M^rj) = (Mw{P) + Y. N &) M f (v)) , (4.4) 

/ 



"Called R2 in the OPP language. 
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Figure 2. Diagrams contributing to the fermionic part of the amplitude (obtained with FeynArts 
[39]). 




y y y y y y 




y y 



Figure 3. Diagrams contributing to the bosonic part of the amplitude (obtained with FeynArts 
[39]). G denotes Goldstone bosons, while u and u ± are the charged ghost and anti-ghost fields, 
respectively. 



where A4w and Mf are scalar form factors of mass dimension -1, and 

= k%k% - (h • k 2 ) g»" . (4.5) 

A4w an d ■M-f are obtained from the diagrams depicted in figures 2 and 3 with the 
Feynman rules of appendix A. Because we work in an arbitrary Rf gauge, also ghost, scalar 
and mixed vectorial-scalar loops contribute to the bosonic part. Notice that diagrams can 
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to 




TO 



s 



Figure 4. Scalar triangle Co(0, 0, mo; m, m, m); dashed lines refer to massless particles, solid lines 
to massive ones. 

be distinguished according to the charge flow in the loop; however, because loops only 
couple to photons, graphs with the same topology but oppositely charged loops equally 
contribute to the amplitude. Schematically, naming the contribution of each diagram after 
the type of particles running in the loop, 



for the bosonic part. Because we work in the i?£-gauge, each bosonic diagram depends 
on the gauge parameter, e.g. Msv — Msv(P,£), but their sum is gauge-invariant, i.e. 
Mw = Mw(P)- 

Even though the final result is finite, UV divergent loop integrals are encountered at 
intermediate steps. In particular, after simplifying reducible numerators, we deal with in- 
tegrals with 2 or 3 internal legs and up to tensorial rank 2. Furthermore, the analytic 
structure of each diagram is characterized by a single threshold, either or m^, or two 
thresholds M w and £,M W . To regularize the infinities, we have used the FDR method 
described in section 2. Instead of calculating the integrals directly, we reduced the ampli- 
tude to scalar integrals by using the FDR version of PV reduction explained in section 3. 
The amplitude, expressed in terms of bubbles, self-energies and triangles (Aq, Bq and Cq, 
respectively), is manifestly gauge-invariant, so that there remain only scalar integrals with 
a single physical threshold to be evaluated. In particular, all divergent integrals cancel 
out, and only scalar triangles must be computed, of the type given in figure 4. Indeed, by 
combining the diagrams, performing the PV reduction, and evaluating the scalar integrals, 
we come by a result consistent with that in [33], which was obtained in DR, i.e. 



Mf = 2M.FFF , 



(4.6) 



as far as the fermionic contribution is concerned, and 



M w = +2 Msss + 2 Mvvv + 2 M S vs + 2 M SS v + 2 M V ss 

+2 Msvv + 2 Mvvs + 2 M V sv + 2 M GGG+ + 2 M GGG 
+2M S v + 2M V s + M S s + Mw , 



(4.7) 



.3 



(An) 2 s w M w . 



2 + 3/3 + 3/3(2 - /3)/(/3) 




Mf{v) 



(Att) 2 s w M w 



2rj i + (i- v )f(rj) 



(4.9) 
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Q + P2 



Figure 5. Vector-scalar-scalar loop diagram contributing to the bosonic part of the amplitude. 



where sw = sin#iy is the sine of the Weinberg mixing angle, and 

arctan 2 ' 
In 



/(*) 



1 



In 2 



2 ( 1+^/1-x+ie 



4 V ~~ l+Vl— x+ie 
is a parametrization of the scalar triangle 
C (0,0,s;M,M,M) = 



l- 



Z 7T 



if X > 1 

if x < 1 , 



(4.10) 



/(*) 



4M" 



(4.11) 



Note that the bosonic form factor in eq. (4.8) contains a constant term, independent of 
the kinematics. FDR automatically leads to its correct value, while in non-gauge-invariant 
frameworks it is necessary to enforce gauge-invariance [25, 26] or momentum routing in- 
variance [23] as an extra constraint to recover it. 

4.1 The W loop contribution 

In this section, we work out in some detail the contribution to the amplitude coming from 
the diagram in figure 5. This stands as an example to illustrate how to use the FDR method 
in practice. In particular, we are interested in showing that, thanks to the global treatment 
of the /i parameter, the integrand within an FDR integral can be algebraically manipulated 
just in the same way as any integrand of a standard 4- or n-dimensional loop integral. 
This fact, together with the FDR shift invariance properties 6 , preserves all cancellations 
required by gauge invariance. 

The analytic contribution of the diagram in figure 5 to the amplitude A4yi/ is given 

by 



Jvl vss 



where, according to the Feynman rules in appendix A, 
e 3 M w {2q + Pl +p 2 ) u (q + 2p2) 



2 sw 



D D 1 D 2 



tip 



:i- e) 

ZMZ) 



(4.12) 



(4.13) 



5 In eq. (4.10), e > is a small imaginary part allowing for the analytic continuation of the result to any 
value of x. 



5 See eq. (2.10). 
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with Ml = d = M 2 , Ml = Mf = £M£ and p n = Ya=i k i- After contracting Lorentz 
indices, q 2 should be promoted to q 2 , as implied by the definition of the FDR integral. The 
usual algebraic manipulations can then be performed on Jvss(q 2 )- Consider, for example, 
the term 



772 



which should equate 



A)£>i£>2 



_1_ | _ dp _ 
D V D 2 DqD 1 D 2 



(4.14) 



(4.15) 



The integrands in eqs. (4.14) and (4.15) give the same result upon FDR integration if they 
generate the same remainder after removing the divergent vacuum configurations. The 
three denominators in eq. (4.14) can be rewritten as in eq. (2.13) 



1 



Do D 1 D 2 



1 



d 2 

— 2 + 



d\ _^ do_ 



<fD 2 q 4 DiD2 q 2 'DoD 1 D 2 



(4.16) 



while only the first term in eq. (4.15) needs an expansion (since the second one is finite): 



1 



D X D 2 



t 



+ 



+ 



tD 2 ' q 2 DiDi 



(4.17) 



The squared brackets indicate the vacuum configurations which need to be subtracted, so 
that the insertion of eqs. (4.16) and (4.17) into eqs. (4.14) and (4.15), respectively, gives 
the expected result 



/ 



Do D 1 D 2 



1 



+ 



do 



D l D 2 DqD x D 2 



lim 



do 



di 



do 



<Fq { —r^ + „— — + — 

I tD 2 q 2 D l D 2 DoD v D 2 



(4.18) 



and one can identify 



1 



D Q D X D 2 



+ 



do 



D X D 2 DoD x D 2 



(4.19) 



at the integrand level. It is important to realize that eq. (4.16) reproduces the last integral 
in eq. (4.18) only if the original q 2 appearing above the three denominators in eq. (4.14) is 
also promoted to q 2 . Otherwise it would give 



lim 

At->0 



& A q 



+ 



q 2 d! 



+ 



q 2 d 



q 2 d2 

tD 2 q i D l D 2 T?D Q D X D 2 



(4.20) 



which differs by a factor m 2 /2 from the correct result. 
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Figure 6. Fermionic loop diagram contributing to the amplitude. 



After simplifying all tensorial integrands appearing in Jyss to irreducible ones, we 
obtain 



"Vss = 4fM (2 + (pi + P2V ^ 2tia ) + 2p% (2 + (p! + »)" C 2thr ) + 



_pj_ 



2 fc r " CfL + (Pi +P2V fe thr - C7f 



.thr 



+ 



2 " BgJ + (ft + ft,)" Khr " ^ 



_thr 



(4.21) 



where the remaining integrals 7 are scalar and tensorial bubbles and triangles in FDR, to be 
further reduced via PV reduction. In eq. (4.21) the subscript 'n thr' indicates the number 
of different thresholds; explicitly 



u 2thr 



7M1 



D (M^)Di(eM^)Z> 2 (CM^) 



■l-'r 



thr 



-°2thr 



R /il.../i r 
"lthr 



7M1 fif^f 



7^1 /yMr 



IT < 



L» (M^)£»!(eM^) 



D (eM^)D 1 (eM^) 

4.2 The fermionic loop contribution 



(4.22) 



The contribution to the decay amplitude mediated by a loop of a charged fermion with 
mass mf is given by 



M F V FF 



2 svv Myy 



Tr [( g + f 2 + m/h^ ^+ ^+ nyh^j + myj] 



7 We put between parentheses the mass appearing in each propagator. 
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with Mo = Mi = M2 = rrif, and where two types of traces containing twice the integration 
momentum appear, namely 



TV[ j/ 7 " il V _] = Tr[(^±/i)7^(^±/i)7^ 



(4.24) 



with d defined in accordance with eq. (2.17). This again allows all the usual algebraic 
manipulations at the level of the FDR integrand; for example 



Tr[ f-yV-f 



[A] 



D Q D 1 D 2 



[A] 



4 q 2 g» v 

1 



4(T ( j[d*q]=-=-+ /[A] 



m 



f 



D D x D 2 



(4.25) 



After simplifying all reducible numerators we obtain 



Mf FF 



S 2 

e rrij r 



s w 



M w \ 



2g» v B + 



+ 4(pi + pa)^ + 4p?C" / + 8C* 4 " } , (4.26) 



where 



■ d 4 g_ ■ ^ ^ 



(4.27) 



After a PV reduction of J3 and the result of eq. (4.9) easily follows. 



5 Conclusions and outlooks 



We made use of the loop-mediated decay H — y 77 to illustrate the key features of FDR, the 
Four Dimensional and free of infinities Regularization/Renormalization approach recently 
introduced in [18]. In particular, we pointed out the mechanisms which lead to an automatic 
preservation of gauge invariance. To this aim, we performed the calculation in a generic 
gauge, showing that, unlike other four-dimensional methods, FDR naturally produces the 
correct rational part of the amplitude, with no need to impose extra constraints. 

The same FDR mechanisms which respect gauge invariance at 1-loop are expected to 
work unchanged with more loops and in the presence of Infrared and/or Collinear diver- 
gences. The explicit verification of these assertions will be the subject of further investiga- 
tions. 
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A Feynman rules 



We draw, in figure 7, the Feynman rules used throughout this work [40]. The tensors V^ up , 
yLLvpcr COU pi j n g constants are given by 

vg p ° = g °f>(p + - P _r + g pp (p+ - qT + g pa (q -p-Y; (A.i) 

y^pc = 2g ^ g ap _ g ^ gP u _ f u gW . (A 2) 

j M w /s w if SVV = HW+W- 
CSVV = \-M W if 577 = ^7 ; ( 3) 



J 2 if 55yy = 77G+G" 

Cssvv - < _ 1/2s ^ . f = ^± 7G t f ; 



1-1 if VS 1 S 2 = ~/G+G- 

Gvs ^ -\ T ^ if VS 1 S 2 = W±G*H ■ (A - 5j 

G denotes Goldstone bosons, while and u^ 1 are the charged ghost and anti-ghost fields, 
respectively. 
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